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We present an analytical method to extract observational predictions about non linear 
evolution of perturbations in a Tolman Universe. We assume no a priori profile for them. 
We solve perturbatively a Hamilton - Jacobi equation for a timelike geodesic and obtain 
the null one as a limiting case in two situations: for an observer located in the center of 
symmetry and for a non-centered one. In the first case we find expressions to evaluate the 
density contrast and the number count and luminosity distance vs redshift relationships 
up to second order in the perturbations. In the second situation we calculate the CMBR 
anisotropics at large angular scales produced by the density contrast and by the asymmetry 
of the observer's location, up to first order in the perturbations. We develop our argument 
in such a way that the formulae are valid for any shape of the primordial spectrum. 

I. INTRODUCTION 

The study of large - scale structure formation represents one of the most exciting research fields in Cos- 
mology. The currently accepted view is that the structures we observe today trace their origin to primeval 
density inhomogeneities, generated by zero point fluctuations in the scalar field responsible for Inflation. The 
evolution of these inhomogeneities, from certain initial values (<< 1) to their present distribution in stars, 
galaxies, clusters of galaxies and so on, involves complex non linear hydrodynamic and gravitational pro- 
cesses. Due to this complexity, it has not been possible to study this evolution from a completely analytical 
perspective §,§, §. 

Given the impossibility at present (and possibly in principle [Q ) to find a general exact analytical solution 
to Einstein equations, our understanding of processes such as the clustering of galaxies is largely dependent on 
numerical techniques. Of the known exact solutions to those equations, most of them presuppose the presence 
of certain symmetries of space-time and/or energy density One of them is the Tolman solution, which 

describes an inhomogeneous Universe filled with pressureless matter and spherically symmetric around a 
point. This solution is completely characterized by two time independent functions / 2 (r) and F(r). This 
last one can be interpreted as the mass contained in a sphere and / 2 (r) as the mechanical energy of a shell, 
both of radius r. 

In spite of its limitations, the Tolman solution has been fruitfully used to study a great variety of effects 
related to the formation of large scale structures: formation of large scale voids, anisotropics in the microwave 
background radiation, possible fractal distribution of galaxies, etc. [Q, ||, [[To), |pT| , |l2| 13 1. But all (or 



almost all) of these studies are based on the numerical integration of the corresponding equations (mainly 
the equation for null geodesies), and often simplified expressions for f 2 (r) and F(r) are employed, sometimes 
with no other reason than to facilitate the numerical calculations. 

The goal of this paper is to present an analytical method to study the non linear evolution of density 
perturbations in a matter dominated Universe, viewed as an instance of the Tolman Universe, laying the 
emphasis in a physically motivated choice for the Tolman functions. Concretely, we shall assume the Tolman 
functions to be such that at early times, in the linear regime, the model be equivalent to a spatially flat 
Friedmann - Robertson - Walker (FRW) matter dominated Universe with growing perturbations. Otherwise, 
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we shall keep the Tolman functions completely general. Observe that in this way the Tolman functions are 
determined by the spectrum of perturbations to the original FRW model. 

We shall focus on simulating in our spherically symmetric model the results of concrete cosmological 
observations such as the redshift - luminosity distance and the number count - redshift relationships. Since 
these observations rely on information carried by light, the analysis centers on the study of the null geodesies 
in the model. We device a perturbative method based on the Hamilton - Jacobi equation for a time like 
geodesic and obtain the null geodesic as a limiting case. 

We apply our scheme to study two situations. In the first one we consider an observer located in the center 
of the Universe and evaluate the number count and luminosity distance vs redshift relationships. For the 
first cosmological test we assume the simplifying hypothesis that the luminosity of galaxies does not evolve, 
which is justified in view of the fact that the redshifts we consider are small (z < 0.08). For larger redshifts 
however, the evolution of luminosity ought to be taken into account Jl4| ). 

In the second situation we consider an observer located away from the center and study the anisotropics 
in the cosmic background radiation due to the presence of perturbations in the photon path, from the last 
scattering surface to the observer's position. 

Throughout the paper a unique matter component of the Universe is assumed. 

We begin our discussion by briefly reviewing the theory of linearized scalar adiabatic perturbations to 
a FRW background, the main features of the Tolman solution to the Einstein equations and the matching 
procedure. In section III we review the cosmological observations N vs z and di vs z and write their form for 
a Tolman Universe. In section IV we begin by sketching our method, we perform the explicit calculations 
for a centered observer and find the general expression for the mentioned observational tests. In section 

V we do the same as in section IV but for a non centered observer and evaluate the anisotropies in the 
CMBR temperature due to the different paths of the photon from the last scattering surface to the observer's 
position. To obtain an idea of the kind of results to be expected from concrete perturbation profiles, in section 

VI we apply our formulae to a particular Tolman Universe, corresponding to a scale invariant spectrum of 
perturbations. Finally we discuss our results in section VII. 

The overall conclusion is that the method we propose affords a simple way to model the growth of nonlinear 
structures such as large scale voids, and may be used profitably in testing competing theories of primordial 
fluctuation generation. 



II. MATCHING A FRW UNIVERSE WITH PERTURBATIONS TO A TOLMAN SOLUTION 



As stated in the Introduction, the accepted framework for structure formation consists in non linear growth 
of perturbations during the matter dominated era of the Universe. The origin of these perturbations are 
the zero point fluctuations of the scalar field responsible for driving Inflation. As the Universe inflates, the 
wavelengths of the modes that build fluctuation grow and become bigger than the Hubble length, leaving 
the region where microphysics can act. At this stage, the amplitude of each mode is frozen and so remains 
until the mode reenters the Hubble sphere. As long as a mode is outside the horizon its evolution can be 
studied within linear theory. We therefore can have an accurate prediction of the amplitudes of all those 
modes whose wavelengths are equal to or bigger than the Hubble length of the considered epoch. Substantial 
growth, moreover, occurs only during the matter dominated epoch of the Universe (modes that reenter the 



horizon during the radiation dominated epoch grow logarithmically at most |15 



In this context, we will characterize the Tolman Universe by matching it to a flat Friedmann Robertson 
Walker Universe with perturbations, at the epoch of equilibrium between the matter and radiation energy 
densities. 

To simplify the matching procedure we only consider modes bigger than (or at most equal to) the Hubble 
size at that epoch. This is certainly justified in hot dark matter models, where smaller fluctuations are 
washed out by Landau damping and free streaming; in more general models a consideration of the transfer 
function would be called for. We shall neglect decaying perturbations to a FRW model, but otherwise do 
not assume any particular spectrum for the perturbations. 

We shall also assume that light traces matter, i.e. a Universe with only one matter component. 
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A. Scalar Perturbations in a Flat FRW Universe 



The most general metric tensor that describes spherically symmetric scalar perturbations in a spatially 
flat FRW Universe is [fil 
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where </>, tf>, B, E are scalar functions of r and £, a prime means partial derivative with respect to r, a dot 
means partial derivative with respect to the cosmic time t and a bar means spatial covariant derivative. 
a = (t/ta) 2/ ' 3 is the scale factor for an unperturbed FRW Universe, where the density is critical, Slo = 1- 

We choose to as the age of the Universe at the time of equilibrium between the matter and radiation 
energy contents of the Universe ((1 + z) eq = 2.32 x 10 4 ), t ~ 5.65 x 10" 4 /i _1 Afpc (the speed of light, 
c = 1, and h is the fudge factor for the Hubble constant, H = lOO/i kms^Mpc^ 1 — » 3 x 10^ i h MpcT 1 in 
natural units). 

The theory of linearized scalar perturbations is treated extensively in Ref. . In the longitudinal gauge 
(B = E = 0), in which the equations take the same form as in the gauge invariant formulation, with a stress 
energy tensor whose spatial part is diagonal (which gives ip = <f>), if we consider adiabatic perturbations with 
wavelengths bigger than the Hubble length, the linearized Einstein equations can be written as a conservation 
law for the quantity 



2^ 5V> 
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where H = a/a. For a matter dominated Universe, the solution is 
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whereby £ = C2(r) follows. The relationship between £ and tp can also be written in terms of the respective 
Fourier modes as 
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Cfc = &Pk/ P represents the departure from homogeneity in the density, caused by that mode of t he p ertur- 
bation, and is directly related to the anisotropies of the microwave background radiation, ST/T ( (llf, p9| ); 
k 3 (k is related to the excess mass SM/M over a region of volume fc~ 3 . From COBE measurements we know 
that k 3 (k — 10~ 5 g(&)j where g(k) is a smooth, order one function of k pfj. 



B. The Tolman Universe 



The Tolman metric for the motion of spherically symmetric dust |2l|] is 



ds 2 = _ R 'l^ dr 2 - R 2 (r, t) (d0 2 + sin 2 9dcp 2 ) (5) 



P{r) 

The Einstein field equations reduce to a single equation pi 
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where R(r,t) > and f 2 (r) and F(r) are two arbitrary functions. F(r) can be interpreted as the mass 
contained in a sphere, and / 2 (r) — 1 as the mechanical energy of a shell, both of radius r. 
If / 2 ( r ) — 1) t ne Tolman equation can be solved explicitly to yield 



R(r,t) 



9 

-F(r) 

4 V ; 



1/3 



(* ~ *o) 2/3 (7) 



Otherwise, we do not have an explicit solution, but the relationship of R to t can still be given in terms of 
a parameter rj as follows: 
/ 2 >1 

R(r,v) = 2[f ^_ 1] (coshr,-l) (8) 



F(r) 

* (r ' ^ = 2[/ 2 (r) - l] 3 / 2 (Shlh 77 ~ ??) (9) 



/ 2 <1 



i?(r ' ??) = 2[/^r) ) -l] (1 ~ COS7?) (10) 

If / 2 — 1 ~ /c7- 2 and F(r) ~ r 3 , the Tolman solution reduces to a dust dominated FRW Universe. In 
this limit, rj becomes the well known "conformal" time. Observe that, in this case f 2 (r) equal to, larger, or 
less than 1, corresponds to flat, open or closed spatial sections respectively. The Tolman solution allows for 
a third arbitrary function of r, namely, an additive constant in the right hand sides of eqs. (Q) and (111]). 
However, we shall make the simplifying assumption of disregarding this extra freedom, which can be shown 
to be of little relevance to the results below. 



C. Matching FRW with Perturbations to Tolman 

Let us now rewrite the perturbed FRW Universe of subsection A in the shape of a Tolman metric. 
To achieve this, we perform an infinitesimal change of coordinates tToiman = ^frw +T(r,t), rxoiman — 
tfrw +K(r,t), and match g™ man = g™ w and Kj° lman = Kfj RW (K tj the extrinsic curvature of a 
spacelike surface) at the surface of matter - radiation equilibrium. Starting from the longitudinal gauge we 
get 

nr,t) = ~^^-C 2 (r)t (12) 

»(r,t) = A^-§^ 8 c/ a (r)* 9 /» (13) 
f(r) = l-frC' 2 (r) (14) 
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We obtain the same expressions for / 2 (r) and F(r) if we start from the synchronous gauge, T and 1Z being 
much simpler. We also note that the expressions for / 2 (r) and F(r) depend only on C^r), the growing mode 
p2fl . In general the r that appears in the FRW expressions is a generic length scale while the one in the 
Tolman expressions is the coordinate distance from the center of the Universe. In our matching procedure, 
the FRW coordinate becomes a radial coordinate. This is not contradictory in view of the fact that FRW 
Universes can be considered a special case of the Tolman solution. 

The relevant growth of perturbations starts at to, the beginning of the matter dominated epoch of the 
Universe. Modes that reentered the horizon before this time at most grow logarithmically, and short wave- 
length modes are suppressed through free streaming and Landau damping | fl9|| . Therefore at to we can 
write the perturbation as a superposition of all modes with wave number smaller than a certain cut-off ko , 
corresponding to the commoving wavelength of the smallest structure that could be formed. We then have 

C 2 (r) = J k 2 dk( k ^ (16) 
o 

where we took the Fourier expansion in spherically symmetric plane waves due to the symmetry of the 
Universe. Cfe is the spectrum of primeval fluctuations, and is to be specified on a t = const, surface. 
Let us write 

C* = jS(fc )e(«, fco) (17) 
where u — k/ko and £ is an order one function. Then C^ir) will take the form 

C 2 (r)~ 7 (fcoKW (18) 



where x = kor and 



{(x) = - I dutt a £(u)=^ (19) 
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The form of 7(&o) depends on the functional form for £(fc), which for the moment we leave completely general. 
In Section VI we will find concrete form for 7(^0) by choosing a particular spectrum for the perturbations. 

^From now on we will employ x, rather than r itself, as radial coordinate. Of course, this coordinate 
change does not affect the form of the Tolman metric. The expressions for F(r) and / 2 (r) become 

n*) = ™[l-57£(*)] (20) 

f( X ) = 1 - y^W (21) 



III. COSMOLOGICAL OBSERVATIONS IN A TOLMAN UNIVERSE 



Having at our disposal an exact solution of Einstein equations, we should in principle be able to follow 
the nonlinear evolution of structures in the Universe in all detail. Our observations, on the other hand, are 
restricted to that part of the Universe that we can "see" , that is, to the past directed light cone with vertex 
at Earth and now. Therefore in order to make contact with directly observable magnitudes, we must first 
discuss how these quantities arc related to the parameters in the Tolman metric, and also how they are 
manifest to us as we look back to the decoupling era from our present location in space time. 

Of the possible observations of cosmological relevance, we have chosen to discuss in detail the dependence of 
luminosity distances and number counts with redshift, and the anisotropy in the temperature of the cosmic 
microwave background. In this section we shall introduce these concepts, discuss the shape of the past 
directed null geodesies in the Tolman Universe, and derive the main formulae to reconstruct the predictions 
of these models concerning redshift surveys and the Hubble law. 
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A. Luminosity distances 



Let us begin our discussion introducing a notion of distance to a cosmological object which can be inferred 
from Earth - based observations. Let dfl be the solid angle subtended by a bundle of null geodesies diverging 
from the observer, and let dS be the cross - sectional area of this bundle at some point. Then the observer 
area distance D a of this point from the observer is defined by |2j| 

dS Q = D 2 Q dn o (22) 

Thus we can find D a if we can measure the solid angle subtended by some object whose cross - sectional 
area can be found from astrophysical considerations. D Q is the same as the corrected luminosity distance 
|p3| and also the same as the angular diameter distance |24|] . For a Friedmann - Robertson - Walker Universe 
the expression for D Q is 

D = a{t s )x (23) 
where and t s indicates time of source emission. For the Tolman Universe we have 

D = R{x,t s ) (24) 



The luminosity distance de is defined as 



where F s is the flux of the source measured in its neighborhood and F is the observed flux. The corrected 
luminosity distance d is defined as 

d = 7rr^ ( 26 ) 

(1 + zj z 

where z is the redshift of the source measured by the observer. By the reciprocity theorem p3[ , we have 
d = D Q , and the luminosity distance di for the Tolman Universe becomes 

dt = R(x,t)(l + z) 2 (27) 



B. Number counts 



One of the cosmological quantities used to study the large scale structures in the Universe is the number 
count of a certain structure (galaxy or cluster) as a function of the redshift z. The starting point to calculate 
this quantity is the number of sources in a section of a bundle of past null geodesies [ p3| , which for a Tolman 
Universe reads [ fic| 

dN = Anv—^—dx (28) 



v being the number density of structures f2]jj 

Fl(x) 



16nM s R/R 2 



(29) 



and where M s is the rest mass of the structure. The number of sources which lie at radial coordinate 
distances less than x as seen by an observer at x = is p0[ 



G 



= *k I ^TW <30) 

LC 

where the integration is made along the light-cone, LC, parametrized by x. For an almost spatially flat 
Universe we can consider that this number (omitting the mass-dependent prefactor) is 

N{x) = F{x) (31) 

This quantity ought to be evaluated along a past null geodesic. 



C. Redshift 



The redshift of a source as measured by an observer is defined in terms of frequencies by 

1 ^emitter /oo\ 

1 + z = (32) 

^observer 

where v e mitter is the intrinsic frequency of the source and ^observer is the value of that frequency detected 
by the observer. If A is an (affine) parameter along the null geodesic between source and observer and t is 
the cosmological time, (|32|) can be written as 

1 | - - (33) 

where x denotes the position of the source and where we take the observer located at x = 0. We shall discuss 
further this formula in next subsection. 



D. Anisotropies in the CMBR temperature 



In the Tolman Universe, when we move the observer away from the center of symmetry, the observations 
that she/he performs will be anisotropic. The physical system on which the most sensitive analysis of 
anisotropies are being performed is the CMBR. In the case of the Tolman Universe there are two sources 
of anisotropies: one related to the fluctuations in the density at the surface of last scattering, AT/T ~ 
(1/3)[AS/S)dec, and the other related to the different paths of the photons from the last scattering surface 
to the non-centered observer, AT/T ~ \z{ff) — z]dec- We will calculate both anisotropies for our model and 
compare one to the other. 



IV. RADIAL GEODESICS IN THE TOLMAN UNIVERSE 



Having derived the generic expressions for luminosity distances and number counts in terms of Tolman 
coordinates and the F and / functions, the only remaining step is to find the evolution of these coordinates 
as we travel towards the past along a null geodesic. To this end, it is convenient to think of a null geodesic 
as the limiting case of a time like geodesic. Time like geodesies have a natural parametrization in terms 
of proper time or some other affine parameter, and the form of the geodesic can be found by integrating a 
Hamiltonian system. We obtain a null geodesic by further imposing the constraint that the Hamiltonian be 
null; the affine parametrization of the time like geodesies then induces a well defined parametrization on the 
null geodesic. 

We can write the action for radial geodesies pa] in Tolman Universe as 



S 



d\{ - 



R' 2 (x,t) 
P{x) 



R 2 {x,t) 



sin 



dX 



(34) 
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where A is an affine parameter for timelike and spacclikc geodesies. The canonically conjugated momenta 
and the Hamiltonian are 

P e =R 2 (x,t)^ (37) 

P^ = R 2 (x,t) sin 2 9^- (38) 
rfA 

H = + R^J) P " + RWrf" + R^{x,t) sin 2 9^ m 
For H — > we have a null geodesic , and in this case the redshift is given by (cfr. eq. ([33|) ) 

(dt/d\)\ x=0 P t \ x =o 1 ' 

We choose the parameter A so that Pt\ x= o = 1, then P t \ x is directly the redshift at that location. 

In order to proceed with the integration of these equations, we should know the expression for R as a 
function of t and x, which entails solving the parametric equations (^) and (^) or ( |l0| ) and [ll]), with known 
F and /. Obviously this cannot be done in closed form unless in the trivial case f = 1. Nevertheless we 
can solve those equations perturbativcly, developing the hyperbolic (or trigonometric) functions around the 
origin (t — ► 0, x ~ 0), the corresponding expressions for R, R~ 2 and R'~ 2 can be read in Appendix A. 

Let us consider the evolution of the perturbations in a neighborhood of the origin. We have to develop 
the functions of x around x ~ 0. For the function F we approximate it by 

4r 3 

F{X) ~ W (41) 
since 7 ~ O(10~ 5 ). From eqs. (^) and (Q) we obtain 

sinh 77 - 77 ~ It l - — ,,V 9 ,J 42 

(of course the same argument applies in the trigonometric case) To achieve a non null value of 77 at the origin 
we have to require that x£'(x) goes to like ax 2 , where a is some constant of order one. In this case eq. 
(H) yields 



t]today 



Hkffit today ) 1/3 ] f^ (43) 



Since 77 decreases as we move down the geodesic, we find that throughout we remain in a neighborhood of 
i] = 0. We are led then to an approximation scheme whereby, in the m — th order approximation, we replace 
the hyperbolic (or trigonometric) functions in the parametric equations by the first nontrivial m terms in 
their Taylor expansion around the origin. For example, at leading order we approximate sinh 77 — 77 ~ 7y 3 / 6; 
at next to leading order we retain also the 77 s term, and so on. In other words, we are neglecting terms 
of order rj 2m+2 /(2m + 2)1 against those of order i] 2m /(2m)\. This approximation is seen to be valid if 
?7 2 <C (2m + l)(2m + 2). rj to day is a measure of the non linear evolution of the perturbation: higher 
Vtoday means higher non linearity and therefore more orders are to be considered in the development of 
the hyperbolic (trigonometric) expressions of t and R. ^From eq. ( fl3| ) we see that increasing r\today means 
increasing kg, i.e. considering a perturbation characterized by a smaller wavelength and that consequently 
is in a more non linear stage of evolution today. 
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A. Leading Order Approximation 



For radial geodesies the Hamiltonian reads 

H = —P? + P 2 (U) 

Starting from eqs. (||) and (^|) and taking the lowest order in the development of the hyperbolic functions 
(77 -C 1) we have 

, 12 r/ 2 (x) - ii 3/2 t , 

, 3 - [/ [ F { X) J (45) 



and therefore 



12 2/3 j pl/3( a .) t 2/3 

R(x,t)~ ^ (46) 



The Hamilton-Jacobi equation for the characteristic function is 

2 -t / Q„A 2 

= —E (47) 



1 fdwv 



with 

F'(x) 



V °W 121/3^2/3(3.)/^) 

and where we used as the principal function S = W + XE. We separate variables by writing W = T(t,p, E) + 
pL(x), with 

(48) 



dT p 2 

*=i E+ W (49) 



As we are interested only in a neighborhood of E = 0, we can expand eq. ( f49| ) 



dT p 1 i 2 / 3 1 9 i 2 

^^^v-s^V"- (50) 



which can be integrated to give 

T-Zpt^z + ^-E- _^_ + ... CD 



•? y-5/3 1 +3 
.1/3 , — E 2 ! 

10 p 24 p 3 



If we define new canonical variables r, p conjugated to -E and p respectively, we have from the characteristic 
function W that 



dW _ 3 t 5 / 3 1 i 3 
~ 10~p~ "12 ? 



r = — = — —-—£;— + ••• (52) 
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from where 



We can write the old momentum P t in terms of the new momenta as 

p _ow _ r~w _ v , v 2/3 i 2 e 

The new Hamiltonian is H = —E. The equations of motion are 

f = dH/dE = -1 
E=p = p = Q 

where a dot means d/dX. As p = const, we can choose it as 



(54) 



(55) 



(56) 



_ ±2/3 
P ~ l today 



so that 



Pt\E=Q 



a + z) 



ttoday 



2/3 



(57) 



therefore when t — t to day we have z = 0. We read the spatial coordinate from eq. ([33]) and as also p = const. 
we take it as 

o - 3t 1/3 



therefore from eq. (&3) we have 



£(i)| B =o = 3^-3tV3 



We can use this last equation to find an expression for F(x). From Eq. (|4S|) and the definition of Vq(x) we 
have 



L(x) 



12 2/3 J-J Tl/S^) 



F 1 / 3 (x)f(x) 
f 2 {x) 



As the second term is 0(7) we can neglect it. So for f(x) ~ 1 we have 



i?1/3 ^ = U2jS L ( X ) = ( 12t to d ay) 1/3 

Finally, recalling Eq. ( [i"l|) the equation for the geodesic becomes 

t i/3 



1 - 



t 1 ' 3 

f l/3 
today _ 



1/3 
today 



(58) 



(59) 



(60) 



where x = 3(fcf|i§t tod ^) 1 / 3 = (!0j«/ 3 ) ^today 

According to eqs. (J5^), (ffq), ( p7[) and (|27| ) the luminosity distance reads 

^ = 3^0^ [(1 + z) - 



(61) 



where we see that for small z, df ~ i/ 1 z, i?o = 2/ 'Sttoday as expected for a flat, matter dominated FRW 
Universe. 
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B. Next to Leading Order Approximation 



Now the Hamilton-Jacobi equation reads 



(dW\ fdW\ 2 1 
\~df) + \~dx~) ?; 2 (x)t 4 / 3 



= -E 



where 



vi(x) 



12 2 / 3 (p(x)-l\ 



2 

Vtoda 



V HI 



at 



2/3 



[FV3( X )]/ V FV3(a;) 
_E and p are no longer constants of motion because in terms of them 

H = -E-p 



today 



2 vi{x) 



10t 2 / 3 



(62) 



(63) 



(64) 



We therefore perform a further canonical transformation to variables E\, pi, pi, ri so that H = —E\. The 
generating function is now 



Wi = Eit + pip - p\Gi(Ei,pi, p,r) 



The old momentum E becomes 



E 



dWi 



so to obtain H — — E\ we demand 



dGi vi(x) 



dr 10t 2 / 3 



We can write formally 



T 

. , _ J_ f , , vi[x(E 1 ,p 1 ,p,T 1 )} 
' 10 7 dT *2/3(£; 1) p 1 ,r') 



(65) 



(66) 



(67) 



(68) 



Under the integration sign we may use the zeroth order relationships among the different variables, namely 
Ei <~ E, pi ~ p 5 and 



dE~dE 



I 



dt'\ E 



£4/3 



dr 



1 t 2 / 3 



2p 



1 + ^ 



(69) 



so 



Gi = 



20pi 



di' 



ttoda 



git 4 / 3 



V! [x(Ei,pi,p,t') 



(70) 



where t — t(r,Ei,pi). Let us now compute the corrections to the several quantities introduced so far. 
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1. Corrections to E 



The value of E valid to this order follows immediately from eqs. 
constraint Ei = 0, leading to 



h and (67) plus the Hamiltonian 



E = — 



Pi vii x ) 

10 t 2 / 3 



where i = t(n, Ei — 0,pi), x — x(ri,pi, E\ = 0,pi) 



2. Corrections to p 



The value of the old momentum p follows from 



P = -tt— = Pi - Pi- 



dp 



Eq. © leads to 



~dp ~ 20 Pl 



ttoda 



dp 



dt' . dx . 



(71) 



where v[(x) — dvi(x)/dx and where the integration is to be performed along a past null geodesic. As v\(x) is 

of first order, dx/dp may be evaluated from the zeroth order, giving dx/dp — %o/3ttod,ay- Therefore, setting 
Ei = we get 



dp 



\E t =0 



•i'O 



20pi 3tl /3 , 

today* 



dt'v'i(x) 



(72) 



and 



P = Pi < 



1 x 
20 -if 1 / 3 

toda yt today 



dt'v'i(x) 



(73) 



3. Corrections to p and t 

The new coordinates are 

dWi n _ „ 2 dG i 



Pi = ^=p-2piGi-pi 1 -^ (74) 
dpi dpi 



dWi 2 dG x 

Tl= m =T - pi dE-i (75) 



For p we have from d70) that at i?i = 



t 

dGi Gi 1 9* 1 [ rU ' > t \ dx \ <-7R\ 

dpi pi 20pi dpi 20pi 7 dpi 
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^From the zeroth order we have p = St 1 ^ + L(x), and consequently dx/dp = and from t = (lOpr/3) 3 / 5 we 
get dt/dp\ T — 3t/5pi. Therefore 

t 

3t If 

p = pl+ 100 t;i(,T) + 20 / dt ' Vl ^ ( 7? ) 

ttoday 

Finally for r we have (cf. Eq. (|75|)) 

t t 
r = r 1 + -^v 1 {x)t I ^ + — ^ / dtt^v[(x)--L f dtt^vxix) (78) 

a ttoday ttoday 



where we have used the zeroth order results dt/dEi\E=o — t 7 ^ 3 /6p 2 and dx/dEi\E=o = x t 5 ^ 3 /Wt]^ ay p 2 . 

4- Redshift, Number Count, Density Contrast and Luminosity Distance 

Having computed the old canonical coordinates in terms of the new ones, it is only a matter of substituting 
them into eqs. ( |54|) and ( p5| ) to obtain the formulae for the number counts (through the auxiliary function 
L(x), recall Eq. (|59|)) and redshifts. For L{x) we get 

5/3 ' 

L(x)=p- 3^/3 + ^_ E t _ =pi + l_J dt ' Vl {x) _ 3t V3 (79) 

o 

Imposing the boundary condition L(0) — 0, we get 

ttoday 

H^)=Ht 1 t l y -t l/3 ) + ^ I dt' Vl (x) (80) 



t 



Substituting the expressions for E and p in Eq. (|5?]) for P t , and imposing the boundary condition P t = 1 
at t = Uoday we get 



2/3 



Now we must find the explicit expressions for dt = dt{x). From the zeroth order we have 

3ttod a y I X \ ^ttoday 



dt = -^^L 1- — \ dx~-^^dx (82) 
x \ x ) xo 

where the last expression is due to the fact that we are interested in a neighborhood of x = 0. Replacing 
this expression and Eq. (p3|) in eqs. (179) and (pl|) we get 



^ \2/3 f 2 

^today \ I 'today ^// / 



p _ / "loaay I -. _ noaay pH / \ 

" ' 1 t / 1 20a ' ' 



j today 



(83) 



(84) 
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and in this case the number count is given by N(x) = F{x) (cfr. Eq. 
by eq. (E9I) omitting the mass factor and in this case we obtain 



AS_ 
T 



2 

^Itoday 

20a 



€ (z)+2 



t 



today 



2/3 



w). The density contrast is given 



(85) 



where 5 = (I2irt 2 ) 1 and where we used S for the energy density in order to avoid confusion with notation. 
Finally, the luminosity distance reads 



df = 3- 



5/3 
today 



t 1 ' 3 

.1/3 
today 



7]today 

20 



107 
3a 



i (x) - 2a* (*) 



t 



'today 



2/3' 



(86) 



where these last three quantities are to be evaluated using x — xq(1 — t 1 / 3 /t]^ ). These expressions are 

completely general, the only assumption on the function £(x) being that x^ (x) goes to zero as x 2 . Therefore, 
given a power spectrum we can find the perturbed background, described by F(x) and f 2 (x) at a given 
initial time to and then evolve it non linearly using the Tolman equations. However, to properly speak about 
non linear evolution, we have to consider higher orders. The procedure to evaluate them is the same as the 
used for the next to leading order. In the next subsection we sketch the evaluation of the second order and 
give the final results. 



C. Second Order Approximation 

In this case the Hamilton - Jacobi equation reads 



fdW\ 2 fdW\ 2 1 



1 



l - Vl {x)t 2 ' 3 + v 2 (x)t i / 3 



10 



-E 



(87) 



where 



v 2 (x) 



12 4/3 



400 [F 1 / 3 (x)}' 



f\x) 1 



i '2 



FV3( a 



12 4/3 



1400 [fV3(ar)]' 



(P(x) 1) 



2 1 



F 



400a 2 i 4 / 3 

today 



11 £" 2 (,)+ :! 



1400a 2 1 4 /3 I x 

today \ 



If we replace in (87) the leading and next to leading order we find again that E\ and pi are no longer 
constants of motion because in terms of them 

H = -E l +ptv 2 (x) 

We perform a further canonical transformation whose generating function is again of the form 

W 2 = E 2 ti +p 2 pi -pjG 2 (pi,Ti,p 2l E 2 ) 
so the first order momentum E\ becomes 

dW 2 2 dG 2 
El --d^- E2 ~ V2 ^ 

To obtain H = —E 2 we proceed as in the previous order and demand 

dG 2 



-v 2 (x) 
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The procedure to calculate the corrections to pi, p\ is the same, although longer, as in the next to 
order case. We therefore show only the final expressions for E\, p\, p\, P t , F 1 ^ 3 (x), AS/S and dg. 

Ei = plv 2 {x) 



V\=Vi\ 



x 



6t 



1/3 



r 

J dtt 2 ' 3 v' 2 {x) 



today 



t 



ttoda 



Pt = 



(^today \ 
~t~) 

3r? 4 



2/3 



1 - 



<2 /3 



2/3 
today 



3^ l^{x) 



800a 2 



A/3 
today 



800a 2 



£ 2 (-) 



2800a 2 

j2/3 



2/3 



t 4 / 3 
A/3 



today today 



t i/3 



4/3 
today 



800a 2 



e 2 (0) 



^/ 3 (*) = (12^)V3 (l - + ^/^H 

k today 



+- 



4000 



IO7 f 

3a 3 7 



da;£ 2 (x) + 



2000 V 3a 3 x 



^Itoday 



* 5/3 , vt 



r(x)-^ + ^-J' 2 (0) + 



2400a 2 s v t 5/3 2400a 2 

today 

4 4 
Vtoday c "2 / \ Vtoday p"2fn\ 

+ 1000a 2 ? 1 J i to<iQ!y 1000a 2 ? W 



A5 = Vfoday 

5 20a 



ttodayj 700a 2 X 2 \ttodayJ 



Vtoday £ (X)(' (X) (^_Y /3 "today ^(x) 

140a 2 a: " UcW 400a2 ^ 2 



\^today J \ttoday J 

Vtoday (X)t" (X) ( t \ 2/3 
200a 2 X \ttodayJ 



2/3' 



I Stodgy t"2 ( x 

+ 400a 2 4 W 



/ t \ 4/3 / t \ 2/3 

\ttoday J \i today J 



+ 



de = 3 



Allay I ^ _t^|_^ + Vtoday ^fW 



today / 



20 V 3a 



2/3 

today , 



t 



4000 V 3a 3 



4/3 

today / 



15 



V today 

1400 



V today 
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£ (*)£ (a:) 13-14 



42/3 44/3 



2800 V 3a 3 



3a 3 



6 + 7- 



.r 



t 2 / 3 
.2/3 

L today 



2/3 1 
today 

t 4/3 



4/3 
today 4 



(95) 



- 9- 



4/3 

today / 



t 



2000 



3a 3 



^Itoday 

' 2400a : 



2^ '(*) — 



5/3 

today 



Vtoday t"2(n\ 

2400a 2 4 1 j 



Jhoday_J>2( s 

1000a 24 W 



today 



^ today p ,f 2(r\\ 

1000a 2 ? 1 j 



and here again, iV(a;) = F(x). Observe that the new corrections are obviously non linear in the perturbations. 

The derivation of these formulae is valid irrespective of the details of the primordial spectrum £(x), and 
so this method may be used to work out in a simple fashion the predictions of several scenarios of primordial 
fluctuation generation. 



V. THE VIEW FROM EARTH 



Due to the spherical symmetry of the Tolman Universe, it is not possible to study anisotropic effects using 
radial geodesies, i.e. an observer located in the center of the Universe, will detect no anisotropies. To study 
the presence of anisotropies in a physical system using the Tolman solution to Einstein equations we have 
to move the observer away from the center of symmetry and solve for it the null geodesies equation. In this 
section we solve the mentioned equation to both leading and next to leading order, and find expressions to 
analyze the anisotropies in the CMBR due to the different redshifts of photons, travelling from the decoupling 
surface to the present. 

Due to the spherical symmetry of the Tolman Universe we do not loose generality if we locate the observer 
at 6 = 0. In this case the Hamilton Jacobi equation for the characteristic function reads 



(dW\ 2 P{x) (dW s 2 



1 fdW\ 2 _ 

R 2 (x,t) \ de 



(96) 



where we used as the principal function S = W 
L(x,p, a) +M(a,9). 



XE. We separate variables by writing W — T(t,p,E) 



A. Leading Order Approximation 



We rewrite Eq. (|96|) as 



fdW\ 



fdW\ 



I fdW\ 2 _ 



\ dt J « 2 i 4 / 3 V dx j vlxH^I* \ 89 j 



E 



(97) 



with Vq = 9£ 2 gJ /a; 2 , and where we have used (|4l]) to evaluate R(x, t) and R (x, t). In this case we have 



dM 
~d9 



voa — » M — vq(x6 



(98) 



vV 



a arccos 



p.r 



(99) 



(100) 
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As we arc interested in a neighborhood of E = 0, we can expand 



dT p 1 t 2 ' 3 1 9 t 2 
dt ~ t 2 / 3 2 P 8 p 3 



which can be integrated to give 



q j-5/3 1 j.3 

T ~Zpt 1 / 3 + -E- -E 2 \ 

1 10 p 24 p 3 



(101) 



(102) 



If we define new canonical variables r, 9 and p conjugated to E, a and p respectively, we have from the 
characteristic function W that 



t _ dW _ 3 t 5/3 1 ^t 3 

T_ aB~io~ _ T2 ^3 



(103) 



p = -7— = — VP 2 a; 2 ~a 2 + it 1 ' 3 - —E^- + ■■■ 
op p 10 p z 



(104) 



dW 
da 



arccos 



px 



^From (104) and recalling that F 1 / 3 {x) ~ (12i tociay ) 1 / 3 x/a;o we get 



FV3(a;) 

(12t today ) 1 /3 P V ^0 



3 +5/3 
10 p^ 



from where we will read the expression for F 1 / 3 {x). 

We can write the old momentum P t in terms of the new variables as 



Pi 



dW 
~dt 



p z p 
^Ta" - ^73 



1 t 2 ' 3 1 2 t 2 

-E E 2 ^r 

2 p 8 p 3 



The new Hamiltonian is H — — E. The equations of motion are 

f = OH/dE = -1 
E=p=p=a=x=0 

where a dot means d/dX. As p = cortsi we can choose it as 



_ j.2/3 

P — l today 



so that 



Pt |b=o= (1 + z) 



j today 

t 



2/3 



(105) 



(106) 



(107) 



(108) 



(109) 



(110) 



therefore when t — t to day we have z — 0. In order to determine a we need to specify the position of the 
observer, xt and the angle of observation, ip. We calculate the observer's position by adjusting the velocity 
of the Earth to the value determined by the dipolar contribution to the background radiation anisotropy, i.e. 



R = rfiday—t' 1 ^ \t=today= 1.666 x 10~ 3 -» x T = 0.833 x 10~ 4 



xo 



(111) 
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The angle subtended by a null ray that reaches xt at t to day is defined through (see Fig. 1) 

P, 



Vq cos(7t — ip) 



p 



from where we obtain 



a = pxT sin tp 



(112) 



(113) 



As p is a constant of the motion, we fix its value by replacing in eq. (104) x — xt, t — ttoday; i- e - we choose 



P — ^ttoday 

Therefore the equation for the geodesies reads 



H cos ip 

x 



X — XQ t 



1 



t l/3 



\ \ t 1/3 



2 1 



today , 



^today J ^ 



XT Xtp 

cos ip H — j 
x 



and the redshift is given by 



(114) 



(115) 



1 + z 



today 



(116) 



B. Next to Leading Order Approximation 



In this case the Hamilton-Jacobi equation for the characteristic function reads 



-E 



where 



Vl 



Wi 



i 4 
10 t 2 / 3 

today 

i 4 

10 f 2 /3 
today 



f\x) 1 



lOat 



2/3 
today 



/ 2 (^) - 1 



lOat 



^ 2 r^j 

2/3 J. 
today 



(cfr. Appendix A) 

Now, £7, p and a are no longer constants of motion because in terms of them 

rr tp 2 V l( X ) « 2 r / \ I \\ 



(117) 



(118) 



(119) 



We therefore perform further canonical transformation to variables E\, pi, oil, p\, \ii T i so that H — —E\. 
The generating function is now 



Wi = E x t +pip + aix -PxSx(p,X, r,p 1 ,a 1 ,E 1 ) - a 2 1 J 7 1 (p,x,T,p 1 ,a 1 ,E 1 ) 
The old momentum E becomes 



E 



dW x 



Ei -pj 



dGi 2 dFi 



dr 



(120) 



(121) 
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so to obtain H = — E\ wc demand 



We can write formally 



dQx _ vi(x) 
Ih ~ t 2 / 3 



dJ- 1 



Qi= dr 



,Vl[ x(p,X,T',Pl,ai, E l)] 

tVZ(E uPU T') 



(122) 
(123) 

(124) 



T\ = / </r 

Ttod 



,wi[x(p,x,t' ,px,ai,E\)] - ui[a;(p,x,r',pi,ai,Ei)] 



^ 2 (/?, X, T r ,pt, ai, Ei) < 2 / 3 (-Bi , pi , r') 
In terms of the zeroth order variables we can write (holding fixed E±, p\) 

1 t 2 l 3 dt 



dT _ d 
dE~ dE 



, p t 4 /3 
dt-^r-Jl + E^- 



't 2 / 3 



ttoday 



dr 



2p 



1 + E 



(4/8 



SO 



2pi 



dt 



l + E^ 



K = 



1 



dt wi[x(p,x,T,pi,ai,£:i)] - ui[x(p,x,r,pi,ai,£'i)] 



(125) 



(126) 



2 Pi J , l\ + e^£- x 2 {p,x,r,pi,ax,E 1 ) 

Ttoday V 

where i = t(Ei,p±,T). Let us now compute the corrections to the several quantities introduced so far 

1. Corrections to E 



(127) 



(128) 



The value of E valid to first order follows immediately from cqs. (121), (122) and (123) plus the Hamiltonian 
constraint E\ = 0, leading to 



2 vi(x) 2 w 1 {x)-v 1 {x) 
where t = t(E± = 0,pi, rll), x — x(E\ = 0,pi, ai, pi, ti). 



(129) 



19 



2. Corrections to p 



The value of the old momentum p follows from 



dW 1 

P = -^r- =Pi - Pi- 



dp 



dp 1 dp 



(130) 



Eqs. (127) and rtl2§) lead to 



ggi _ 1 
dp 2pi 



dt 



ttoday V P 



i . .dx 

= v x\ x )-r 
t'W dp 



(131) 



flFi 1 
dp 2pi 



wi(x) - ui(x) 



dx 
r/p 



(132) 



Setting Ei = 0, we have from the zeroth order 



. - 3* 1 / 3 p - St 1 / 3 _ y/p*x 2 - a 2 



dx 



dp , y(^ + - 



v Pix 



(133) 



and consequently 



p = Pi 



Pi 



dt 



I ; = ,; ( x) ^ 2 - Q " 

today V P l 



(134) 



t /4/3 



ttoday V Pi 



wi(x) - vi(x) 



^p\x 2 - a{ 

VqPiX 



3. Corrections to a 



The value of the old momentum a follows from 



dW x 
a = — — = ai 
OX 



(135) 



a result that could be guessed due to the fact that the perturbations do not depend on the angles. 



4- Corrections to p and t 



The new coordinates are 



9Wi 2 dQ x 2 dT Y 

Pi = -r. — = P~ 2piGi -Px^ aj- — 

dpi dpi dpi 



(136) 



At Ex = we have from eqs. (127) and (128) that 
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dQi 1 f , , , 1 , . dt 1 [ „> , , dx 



, dtvi(x) + — vi(x)— + — / dtv^x) — (137) 
dpi 2p\ J 2pi dpi 2pi J dpi 

Ttoday 1 today 



dTi J_ 

dpi 2p\ 

1 

+ 2p~i 



dt' 



, ( wi {x) — vi (x) \ 1 f wi{x) — vi(x) \ <9i 



2pi 



9p] 



+ 



(138) 



ft oday 



, d ( u>i (x) — vi (x) \ dx 



dx 



dpi 



T~t oday 



^From the zeroth order, if E = 0, we have that t — (10pr/3) 3 / 5 , so that dt/dp = 3t/5pi. Also from this 
order 



Therefore we have 



(p-3tV 3 ) ^ «2 ^ dx 



p V — 53 — + F 



(139) 



p = pi + - / dtvi(x) + — vi(x)t + 



3 a 2 f wi(x) — vi(x) 

To pf 



<t oday 

t 



^■toda 



wi(x) - vi(x) _ a\_ 
x 2 2p\ 



dt 



v^x) 



't oday 



(140) 



2p\ 



dt 



w^x) - v^x) \ a\ f wi(x)-vi(x) 



-T I dt- 

Pi 



T today 



In orde r to compute the integrals, it is convenient to parametrize the null geodesic with the coordinate x, 
from eq. (115) we get 



dt 



1 



3pttoday 

Therefore the corrections computed above read 



PiXj. — d\ — Jp\x 2 — a\ 



2 .2/3 

VottodayPlX 

\/p\x 2 — Ct\ 



dx 



P 



, Pi ,2/3 , x . .2/3 (Wi(x) ~Vi(x)) 
Pi + ^todayMx) + ^t tQday - 



(141) 



(142) 



,2/3 f, x . , 3t , . 3a? wi(x)-Vi(x) 
P = Pi~ t t L v PiVo J dx __^ (a; ) + - Vl (x) + ^t - + (143) 

X X 

. a\ ,2/3 f dx , "i, 2/3 f dx ( w t (x) - Vi (x) \ 

7pW=^ Mx) + ^ W'°/ ^2,2 _ a2 { x ) + 

Qq_ 2/3 f dx w[(x) - v[(x) of 2/3 f dx wi{x)-vi(x) 
2p3 ho day V0j _ ^ x2 p f todays J _ ^ 
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where x is to be read from the leading order and where v%(x) and w\{x) are given by eq. ( |118|) . 

Now we are ready to evaluate the next t o le ading orde r co rrections to the redshift and F 1/i (x). For the 
redshift we have to replace eqs. ( |129| ) and (142) into eq. (107), keeping only the leading and next to leading 
order terms. We get 



Pi 

£2/3 



1 - 



Vtoday 

20a 2 



t 2 ' 3 \ 

2/3 I 

today / 



2 

^toda 



20a 



y 9 . 9 



t 2 ' 3 

.2/3 
today / 



(144) 



where we have used a\ = p\Xt sin tp. 

In order to evaluate F 1 ^ 3 {x) we recall expression (|106|) and replace in it eqs. (143) and (129) obtaining 



F" 3 (i) = (12( M „»)" 3 



u 



H cos tp t-7- h I(x) 



,T 



1/3 
today 



— sin 



(145) 



where 



/(as) 



Vtoday /IO7 



20 



3 a 



x 2 — x|n sin 2 ip 



-4 (x) 



Vtoday /IO7 2 . 2 



rf.E 



x 2 — x|n sin 2 y 



Vtoday /IO7 2 -2 



x 2 — iy sin 2 y 



Vtoday /IO7 4 . 4 



e (*)- 
e' g>) 



x 2 — sin 2 y 



(146) 



(147) 



In both equations (144) and (145) we have imposed as boundary conditions P t — 1 and F 1 / 3 {x) 

(l2t to day) 1/3 X T /x at X = X T , t = t to day 

As earlier, our formulae are valid for any choice of £(x). 



VI. AN EXAMPLE: SCALE - INVARIANT PRIMORDIAL PERTURBATIONS 

In order to apply our developments to a concrete case, we must now specified Cfc m Eq. (|l6|). In the 
absence of a specific model for the spectrum of a primeval fluctuations, it is often assumed that, at horizon 
crossing = pk^ n ~ 4 ^ 2 . On at — const, surface the corresponding expression is 

|a I =0 n ~ 3 V 2 (148) 

where n is known as the spectral index. This index is determined by the equation of state ui = p/ p at 50 — 60 
e-folds before the end of Inflation, and its expression for scalar perturbations is n = 1 — 3(1 + u>) + dln(l + 
u))/dlnk, where k is the wavenumber p6| |. Strict exponential (de Sitter) expansion corresponds to n = 1, 
the so-called Harrison - Zel'dovich spectrum. But in any realistic inflationary model, the expansion rate 
must slow down near the end of Inflation in order to return to FRW expansion, and in this case n =/= 1 is 
expected |26| . The specification of the spectrum is concluded by also fixing the phases of the amplitudes for 
the different modes, which in inflationary models are usually random Gaussian variables; for simplicity, we 
shall disregard this, assuming |^| = Cfe instead. Of course, applying our formulae to the more general case 
of nontrivial phases involves no difficulties of principle. 
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With the form (U4q), we have 



and 



7 = 



(n+3)/2 



(149) 



m = l - f 

X J u 





(n + l)x 



(n+l)/2 



sin ux 



(150) 



1 n 



-ix 



where iiq 2±£; 



% From eqs. @ and (§T|) we 
The effect of the perturbations 



±iar) are the degenerate hypergeometric function |2 
see that in general F(x) and f 2 (x) will be oscillatory functions of x. 
at the beginning of the matter dominated epoch is a mass redistribution, and the creation of negative 
curvature regions (f 2 (x) > 1), that will expand forever, and of positive curvature ones {f 2 {x) < 1), that 
will eventually collapse. This distribution of matter will evolve non linearly from its initial state at to, its 
evolution being described by the Tolman equations. We normalize the spectrum by considering the amplitude 
of the anisotropy of the cosmic background radiation corresponding to a scale of the size of the horizon at 
decoupling (X d ~ 0.13 Mpc) |§. 

(n+3)/2 



kd(k d — lnk d 



(n+3)/2 _ 



7 



6.6 x 10" 



(151) 



where we used k^Ckd — 3AT/T. By virtue of the linear relationship between the angles that the scales 
subtend today and their sizes, we have kg/ho ~ 0a/0d, where 9d is the angle subtended today by a scale of 
the size of the horizon at decoupling, 8^ ~ 1.4° — 2°, and 8 is the angle subtended today by the wavelength 
corresponding to fc , 8 ~ 45' — 1°. For concreteness we choose 8o ~ 48', 8^ ~ 1.4°. 

With the form (152) for £(x) we have a = 1/12. If we take fco ~ 918/iMpc -1 (which corresponds to a 
wavelength about six times bigger than the horizon size at equilibrium), to ~ 5.65 x 10 _4 /i _1 M|>e, t to day — 
lQ§%hr l Mpc and 7 ~ 2.02125 x 10" 4 , we get r] ~ 1.78 (cfr. eq. @). 

The model becomes particularly simple when n = 1, whereby eq. (|150|) becomes 



(152) 



This function is to be used in expressions found in the previous sections. 



A. Density Contrast, Number Counts and Luminosity Distance 

In Fig. 2 we have plotted the density contrast AS/ 6 vs z for the computed second order. As the matching 
conditions tell us nothing about the sign of f 2 (x) — 1, we chose it so that the perturbation evolved to form a 
void. We see that the profile is of non compensated kind. After the non linear evolution, the density contrast 
in the wall is AS/8 ~ 0.06 and in the center of the void AS/6 ~ —0.44. If we define the radius of the void 
as the distance between its center and the point where A6/6 = (which corresponds to z ~ 0.024), we get 
R ~ 72h~ 1 Mpc, taking into account the results of Fig. 4. 

In Fig. 3 we have plotted lnJV(x) vs lnz. Dotted line corresponds to an unperturbed Universe, full line 
to one with second order perturbation. We see that the number count increases with z as we move over the 
increasing density zone. For z ~ 0.05 the number count begins to equal that for a FRW Universe. 

In order to analyze the departures from the linear Hubble law due to the presence of the perturbations, 
we have plotted in Fig. 4 the z vs dg relationship, for a flat FRW Universe (dotted line) and for Tolman 
Universe (full line). Due to the small redshifts considered, the slope of the curves can be considered a 
measure of the Hubble constant, which in our units is Ho — 2 (the true value of the Hubble constant being 
Ho = 2/3ttoday)- We see that the slope of the curve corresponding to the perturbed Universe is smaller 
than the one corresponding to FRW, up to z ~ 0.013. From there on the slope starts to increase to finally 
coincide with the unperturbed one as expected. 
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B. Anisotropies in the CMBR 



In our Tolman Universe, we find two sources of CMBR anisotropies, one of which is due to the fluctuations 



in the density, AT/T oc AS/S (AS/ 6 given by eq.(|S5|) evaluated using eq. (115)). The origin of the other is 
that being the Earth away from the center of symmetry, the light that was emitted from different points of 
last scattering surface will reach the observer with different redshifts, depending on the angular position of 
the source point, i.e. we have a contribution to the CMBR anisotropies of the form z(9) — C, where Q is the 
mean redshift of the decoupling surface. We find that in our case the last contribution outweights the first 
one by a factor of 10 7 : we calculated AS/S and z(9) — Q and after substracting the dipolar contribution due to 
the movement of the Earth (AT/T w 1CU 3 ) we were left with a quadrupolar contribution of AT/T w 10 -13 
and AT/T w 10 -6 , respectively, for big angular scales. The fact that the contribution from the density 
fluctuations is so small is not surprising in view of the small cut-off kg . 

VII. DISCUSSION 

In this paper we have developed an analytical method to study the non linear evolution of adiabatic 
perturbations in a matter dominated Universe. We build the initial profile of the fluctuation on a t — const. 
surface (that corresponds to the beginning of the matter dominated era, when structure formation begins) 
as a superposition of all those modes whose wavenumber is smaller than a certain cut-off fco- This cut-off 



may correspond to the mode that survived the Meszaros effect 1 15 , Landau damping or free streaming J_9| . 

^From this point on we can follow the non linear evolution of the perturbations by means of the exact 
solution to Einstein equation for a pressureless, spherically symmetric Universe, namely the Tolman Universe. 
Since this solution is isotropic with respect to one point, the Tolman solution can be applied only to analyze 
the formation of a structure (for example a void or cluster) with the proper symmetry. This restriction in the 
applicability of the Tolman solution to the study of structure evolution is compensated by the advantages of 
having an exact solution at our disposal. 

Formally, the basis of all calculations for observable quantities lays in the resolution of the equation for 
the null geodesies, which in general is not trivial and is carried out by numerical methods. In this work we 
presented a perturbative method to do this calculation analytically. We put our calculations in a physical 
framework by the way in which we build the perturbations. We carried out the calculations only to second 
order, which means a mild non linear evolution, but the extension to higher orders is straightforward. 

In the last section we obtained a glimpse of the results to be obtained by assuming a scale free spectrum 
of initial perturbations. In view of the simplicity of the calculations involved, we find that our results give a 
remarkable approximation to the structure of some of the largest known voids, such as Botes'. Quantitatively, 
our results are accurate only up to an order of magnitude (in fact, the Botes void is surrounded by a wall 
whose density contrast is AS/S ~ 4, and its radius is of about 30Mpc p9[). The difference, of course, could 
be reduced by a better choice of the cutoff wavelength and, most importantly, by carrying the computation 
to higher orders. 

With regard to future work, the most important feature of the method we propose, besides its analytical 
and relatively simple character, is that it may be applied for any form of the primordial spectrum. Thus it 
becomes only a matter of plugging in one's favorite theory of primordial fluctuation generation (as reflected 
in the particular form of ^ (Eq. (|16|)) to easily obtain (rough) testable predictions of that model. While we 
have used here a (unrealistic) scale invariant spectrum for demonstration purposes, the calculation is equally 
simple with red, blue, or more sophisticated alternatives. 

We therefore believe the methods we describe in this paper will be an useful tool in the delicate task of 
sorting between the manyfold fluctuation generation scenarios now available. 
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IX. APPENDIX A: EXPANSION OF THE RADIAL TOLMAN FUNCTION 
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X. APPENDIX B: EVALUATION OF Q 

We evaluate fi today as follows. In a neighborhood of the origin, the evolution of the Universe is nearly 
like a FRW one, therefore 

n = l+ H^)= l + lk) (156) 

where H is the Hubble constant, a(t) the expansion factor of the Universe and JC the curvature. We have to 
extract the expressions for a(t) and K. from the Tolman equations (|l4|), (||) and (||). We define 

/C = -2^^ I *=° (157) 

from where we get K = (10/3)7^ (0) ~ —5 x 10~ 5 . To obtain d(t) we calculate R — (dR/dr])(dr)/dt) from 
Eqs. (||) and obtaining 



*(.,„- vT^Fi^-^f,-^ (158) 



where a = |£ (0)|. We therefore have 



;10 sinh?7 . . 

^ = \l^ a ^TT (159) 



Replacing Eq. (159) in Eq. fll56|) we get 



25 



2 

Evaluating rj from Eq. with t = t to day we get fttoday — 0.49 at x = 



= 1- tanh 2 (2) (160) 
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FIGURE CAPTIONS 

Fig. 1: Angle subtended by the null geodesies that arrive at the observer's position xt- 

Fig. 2: Density contrast AS/ 5 vs redshift z, to second order in the perturbative development. The profile 
corresponds to a non compensated void. 

Fig. 3: In AT vs In z, to second order in the perturbative development. Dashed line corresponds to flat 
FRW Universe, full line to Tolman Universe 

Fig. 4: Redshift z vs luminosity distance dl, to second order in the perturbative development. Due 
to the small redshifts considered, the slope of the curve is a measure of the Hubble constant. Dashed line 
corresponds to flat FRW Universe, full line to Tolman Universe. 
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Fig. l 



Fig. 2: Density Contrast Vs Redshift 




